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Many features of spin models can be interpreted in geometrical terms by means of the properties of well defined 
clusters of spins. In case of spontaneous symmetry breaking, the phase transition of models like the q-state 
Potts model, 0(n), etc., can be equivalently described as a percolation transition of clusters. We study here 
the behaviour of such clusters when the presence of an external field H breaks explicitly the global symmetry of 
the Hamiltonian of the theory. We find that these clusters have still some interesting relationships with thermal 
features of the model. 



1. INTRODUCTION 

Since the early 40's jij the interest of several 
physicists studying critical phenomena has fo- 
cused on the investigation of geometrical struc- 
tures of the system whose interplay could repre- 
sent the mechanism for the occurrence of a phase 
transition. 

A breakthrough in this course of studies was 
the discovery of Fortuin and Kasteleyn (FK) [|| 
that the g-state Potts model with vanishing ex- 
ternal field can be reformulated as a geometrical 
model, with the spin configurations turning into 
clusters configurations by connecting nearest- 
neighbouring spins of the same orientation with 
a bond probability pb = 1 — exp(— J/kT) (J is 
the Potts coupling, T the temperature). A re- 
markable consequence of this mapping is that the 
FK clusters percolate at the critical point of the 
thermal transition and, in case of second order 
phase transitions, the percolation exponents co- 
incide with the thermal exponents. This result 
has been since then extended to a wide variety 
of models, including continuous spin models like 
0(n) 1. 

The absence of an external field endows the 
Hamiltonian of all these models with a global 
symmetry, which is spontaneously broken by the 
state of the system at low temperatures. The 
spontaneous breaking of such global symmetry is, 
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ultimately, the reason of the phase transition. If 
we switch on an external field H, the Hamilto- 
nian of the system breaks explicitly the global 
symmetry one has for H — 0, and this has deep 
consequences on the possibility for the model to 
show critical behaviour. One has, in general, two 
possible cases, according to the order of the phase 
transition when there is no external field: 

• The model undergoes a second order phase 
transition for H = 0. 

• The model undergoes a first order phase 
transition for H = 0. 

Here we want to see whether the clusters that 
describe the critical behaviour at H = have 
any relationship with thermal properties of the 
model also in presence of an external field. For 
this purpose we will analyse separately the two 
cases listed above. 

2. SECOND ORDER PHASE TRANSI- 
TION FOR H=0 

If the system undergoes a second order phase 
transition for H = 0, for any Hy^O the partition 
function of the system is not singular and one has 
at most a rapid crossover from the ordered to the 
disordered phase. 

From the renormalization group ansatz of the 
magnetization as a function of h = H/kT and 
t = (T — T c )/T c one derives simple scaling laws 
for t, /i<g;l. Particularly interesting is the so called 
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pseudocritical line, which is given by the temper- 
ature t x at which the susceptibility \ peaks for a 
given h. The equation of the pseudocritical line 
is given by 

t x oc h 1 "", (1) 
where j3 and 6 are critical exponents. 
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Figure 1. Comparison of the Kertesz and the 
pseudocritical line for the 2D Ising model. 

On the other hand, the FK clusters^] perco- 
late for any H^O: the corresponding geometrical 
variables show divergences as in the case H = 
and percolation remains a genuine critical phe- 
nomenon. The line of percolation thresholds is 
called Kertesz line 

We want to compare the Kertesz and the pseu- 
docritical line. Fig. [I] shows the two lines for the 
2D Ising model. In the double logarithmic scale 
of the plot, they look parallel to each other. A fit 
of the four points of the percolation line leads in- 
deed to the same scaling law as Eq. |l|; the relative 
exponent is n = 0.534(3), in excellent agreement 
with 1/(35 = 8/15 = 0.533. We have also found 
that the same result holds for the 3D 0(2) model 
H , so that it is likely to be general. 



3. FIRST ORDER 
TION FOR H=0 



PHASE TRANSI- 



In the case of a first order transition for H = 0, 
the model shows a discontinuous phase change 

2 The result is valid for any model which can be mapped 
onto a cluster model, e.g. O(n). 



also in the presence of an external field, as long as 
H is smaller than some critical H c . One has then 
a line of first order phase transition thresholds, 
from H = 0, T = T c (0) to H = H c and T = 
T C (H C ) (endpoint), where the transition becomes 
continuous and the exponents are conjectured to 
be the Ising exponents. 

We have examined the case of the 3D 3-state 
Potts model ||, which is subject of intense inves- 
tigations because its phase transition is closely 
related to the dcconfincmcnt transition of finite 
temperature QCD |t]]. The critical value of the 
field H c was recently determined with great pre- 
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Figure 2. Time history of M and P at a point of 
the first order transitions line. The lattice size is 
70 3 . 



In order to see how the FK clusters behave 
for H^0 we have analyzed the time history of 
the percolation order parameter, the percolation 
strength P, which is the probability that a ran- 
domly chosen site of the lattice belongs to a per- 
colating cluster. Fig. || shows the time history 
of P compared to the one of the magnetization 
M at a point of the first order phase transitions 
line, quite close to the endpoint. From the fig- 
ure we can see that P makes a jump from zero 
to a non-zero value: there is then a discontinuous 
geometrical transition between a percolation and 
a non-percolation phase. We repeated the proce- 
dure for several points of the first order transi- 
tions line, obtaining any time the same result. 
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For H = H c we expect that the percolation 
transition of the FK clusters becomes continu- 
ous like the thermal one. It is then interesting to 
check whether the two thresholds coincide and, in 
this case, whether the percolation exponents be- 
long to the 3D Ising universality class as the ther- 
mal exponents. Fig. || shows the percolation cu- 
mulant, that is the fraction of configurations with 
a percolating cluster, as a function of f3 = J/kT 
for different lattice sizes. The crossing point of 
the curves is the critical point of the geometri- 
cal transition and it is in good agreement with 
the thermal threshold determined in || (vertical 
dashed lines in the figure). 
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Figure 3. Percolation cumulant at H = H c for FK 
clusters as a function of /3 for four lattice sizes. 

The height of the crossing point is a universal 
number: the figure shows that it considerably dif- 
fers from the values corresponding to the 3D Ising 
and random percolation (RP) universality classes 
(horizontal lines in the figure). A finite size scal- 
ing analysis of the cluster variables at the critical 
point leads to the following values for the perco- 
lation exponents: /3/v = 0.32(3), 7/1/ = 2.32(2), 
v = 0.45(3). As expected, these exponents belong 
neither to the 3D Ising nor to the RP universality 
class. 

4. CONCLUSIONS 

We have seen that the clusters which describe 
the critical behaviour of thermal models in ab- 
sence of an external field H maintain a close rela- 



tionship with thermal features of the models also 
when Hj^O. 

If the spontaneous symmetry breaking transi- 
tion is second order, in the limit of small fields the 
Kertesz line is described by the same function as 
the pseudocritical line of the model, a power law 
with exponent 1/(36. 

If the spontaneous symmetry breaking transi- 
tion is first order, the line of thermal first order 
phase transitions is also a line of first order per- 
colation transitions for the clusters. In contrast 
to the magnetization, the percolation strenght P 
is a good order parameter for the geometrical 
transition all along the line. The endpoint is a 
percolation point, but the corresponding critical 
exponents are neither in the Ising nor in the ran- 
dom percolation universality class. That suggests 
that one might need to change the cluster defini- 
tion, eventually introducing an explicit field de- 
pendence in the bond probability ps, so to obtain 
the correct critical exponents of the thermal tran- 
sition at the endpoint. Such definition could turn 
out to be fruitful also in the crossover region. 
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